Dynamical constraints on phase transitions 



Klaus Morawetz 1,2 
1 LPC-ISMRA, Bid Marechal Juin, 14050 Caen, France 
2 GANIL, Bid Becquerel, 14076 Caen Cedex 5, France 

The numerical solutions of nonlocal and local Boltzmann kinetic equations for the simulation of 
central heavy ion reactions are parameterized in terms of time dependent thermodynamical variables 
in the Fermi liquid sense. This allows one to discuss dynamical trajectories in phase space. The 
nonequilibrium state is characterized by non-isobaric, non-isochoric etc. conditions, shortly called 
iso-nothing conditions. Therefore a combination of thermodynamical observables is constructed 
which allows one to locate instabilities and points of possible phase transition in a dynamical sense. 
We find two different mechanisms of instability, a short time surface - dominated instability and 
later a spinodal - dominated volume instability. The latter one occurs only if the incident energies 
do not exceed significantly the Fermi energy and might be attributed to spinodal decomposition. In 
contrast the fast surface explosion occurs far outside the spinodal region and pertains also in the 
cases where the system develops too fast to suffer a spinodal decomposition and where the system 
approaches equilibrium outside the spinodal region. 



I. INTRODUCTION 

The collisions of heavy ions around the Fermi energy 
and the description of multifragmentation phenomena 
fills an enormous literature. Mostly the multifragmen- 
tation is attributed to a hypothetical liquid - gas phase 
transition which is partially supported by mean field con- 
siderations in equilibrium where the nonlinear density de- 
pendence of the interaction energy leads to a liquid - gas 
like first order phase transition. Therefore, the phenom- 
ena has been investigated in terms of a spinodal decom- 
position. However this straightforward picture is over- 
shadowed by at least two serious drawbacks. First we 
have to deal with finite systems, where the phase transi- 
tion appears modified and less pronounced than in infi- 
nite bulk matter. And secondly, we have to face the fact 
that the process evolves under extreme nonequilibrium 
conditions. For a critical discussion of recent models on 
multifragmentation see [H. We want to investigate the 
latter two points here and will use a microscopic approach 
which allows one to describe the time evolution of the 
one-particle distribution function including binary corre- 
lations. We will suggest a possibility of analyzing phase 
transitions in terms of time dependent thermodynamical 
variables and will be able, in this way, to see signals of 
instability in nonequilibrium and finite systems. 

The standard treatment to investigate basic features 
of multifragmentation processes is performed in terms 
of fluctuation analysis starting from the Landau equa- 
tion |-|] or BUU equations |,§. Obser ving that these 
kinetic equations do not lead to enough fluctuations to 
describe multifragmentation, additional stochasticity has 
been assumed and incorporated resulting in Boltzmann- 
Langevin pictures [fz| — p~2 1 . The large time scale of fluc- 
tuations has been analyzed in ||l3| . It is found that the 
large time evolution of the system is guided by coopera- 



tive effects and fluctuations in a universal manner. The 
crucial role of collision rate has been pointed out in that 
it enforces the diffusive regime. 

We will adopt here a straight microscopic picture of ki- 
netic theory without additional fluctuations of Langevin 
sources. While this is perfectly microscopical controlled 
we have to leave out the possibility of describing fragment 
production. In contrast we will investigate the thermo- 
dynamical trajectories as arising straight from the solved 
kinetic equation as a Fermi liquid. Therefore no coa- 
lescence or other cluster creation mechanisms are used. 
This allows to restrict to the single particle distribution if 
the two - particle correlations are included in the collision 
integral. This is performed in the frame of the nonlocal 
kinetic theory. Since we want to study the dynamical 
constraints of phase transitions as necessary but not suf- 
ficient conditions we can expect already from the kinetic 
theory an answer whether the system will undergo spin- 
odal decomposition or other forms of decomposition. In 
fact we will demonstrate that there is a dominant surface 
emission at higher energies than the Fermi energy while 
the spinodal decomposition can be accessed only for ener- 
gies lower or equal the Fermi energy. For higher energies 
the system evolves too fast through the spinodal region 
to be influenced sufficiently by spinodal decomposition. 

There are two experimental hints for two different 
regimes of instability in heavy ion collisions around Fermi 
energy. The first one concerns the emulsion data recorded 
in the experiment by F. Schussler et. al. which has 
been considered in |14| . There the fragments with charge 
Z > 2 has been grouped into two different velocities, one 
around 0.16c and the other with 0.25c. A possible in- 
terpretation has been advocated that the higher velocity 
group comes from fragments emitted at an early stage 
from the surface. The TDHF calculations seemed to sup- 
port this picture. 
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A second experimental signal comes from the produc- 
tion of hard photons as measured by the TAPS collabo- 
ration The extracted source sizes by HBT interfer- 
ometry has found to be too large if not two sources are 
assumed. Moreover the calculated photon spectra shows 
a clear prompt source of hard photons besides later ther- 
mal photons. The latter second source vanishes for inci- 
dent energies larger than 60MeV. This indicates already 
that there is a transition between two mechanisms of par- 
ticle production and instability if the bombarding energy 
exceeds 50-60 MeV. 

We will show that indeed there can be identified two 
mechanisms; at short times a surface dominated emission 
and at later times a volume dominated spinodal decom- 
position. For low energies we will find that the volume 
spinodal effects are visible while for higher energies only 
the surface emission survives. 



II. KINETIC DESCRIPTION 

We use for the description the recently derived nonlocal 
kinetic equation Jl6[ ] for the one - particle distribution 
function 
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with Enskog-type shifts of arguments Pq| : ft = 
f a (k,r,t), h = A(p,r-A 2 ,t), h = /a(k-q-A K ,r- 
A 3 ,i-A t ) , and / 4 = / h (p+q-A A - ! r-A 4 ,i-A t ) . The 
effective scattering measure, the T-matrix is centered in 
all shifts. The quasiparticle energy e contains the mean 
field as well as the correlated self energy. 

In agreement with Jl^Jl^], all gradient corrections 
are given by derivatives of the scattering phase shift 
= ImlnTj(ft,k,p,q,i,r), 
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After derivatives, A's are evaluated at the energy shell 
ft — > £3+64. Neglecting these shifts the usual BUU 
scenario appears. 

The A's in the arguments of the distribution functions 
in (Q) remind the non-instant and non-local corrections in 



the scattering-in integral for classical particles. The dis- 
placements of the asymptotic states are given by A2.34. 
The time delay enters in an equal way the asymptotic 
states 3 and 4. The momentum gain Ak also appears 
only in states 3 and 4. Finally, there is the energy gain 
which is discussed in p9[ . These nonlocal corrections to 
the usual Boltzmann equation are a compact form of gra- 
dient corrections. It ensures that the conservation laws 
contain besides the mean-field correlations also the two 
particle correlations. 

Despite its complicated form it is possible to solve 
this kinetic equation with standard Boltzmann numer- 
ical codes and to implement the shifts |2(| . Therefore we 
have calculated the shifts for different realistic nuclear 
potentials [^l). The numerical solution of the nonlocal 
kinetic equation has shown an observable effect in the 
dynamical particle spectra of around 10%. The high en- 
ergetic tails of the spectrum are enhanced due to more 
energetic two-particle collisions in the early phase of nu- 
clear collision. Therefore the nonlocal corrections lead to 
an enhanced production of preequilibrium high energetic 
particles. 

Besides the nonlocal shifts and cross section which has 
been calculated from realistic potentials we adopt here 
the view that the selfenergy e is parameterized in terms 
of the Skyrme potential for which we use a soft potential 
of the form 
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For a derivation of collision integrals and the Skyrme 
potential from the same microscopic footing, see [^2) . 



A. Balance equations 

By multiplying the kinetic equation with l,p, e one 
obtains the balance for the particle density n, the mo- 
mentum density J and the energy density £ . Without 
nonlocal corrections the collision integrals vanish for the 
density and momentum balance and we get the standard 
balance equations for the quasiparticle parts 
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with the quasiparticle density, the current and the mo- 
mentum tensor 
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where the quasiparticle energy is given by 
dp p 2 
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and the pressure is as usual 

P = -Vu. 
3 

The quasi particle energy of the system varies as 
dp S£ qp 
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and since we adopt the parameterization of quasiparti- 
cle energy (^), the quasiparticle part of the total energy 
density reads 
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Please note that besides the mean field (||) we have also 
a Born correlation term £ Born coming from the second 
term of (0), see p3|, 
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The balance of the quasiparticle part of the energy den- 
sity reads from the kinetic equation 



<9£qP(r,t) d_ 
dt dr 



The correlational parts of the density, pressure and 
energy are coming from genuine two-particle correlations 
beyond Born approximation which are also derived from 
the balance equations of nonlocal kinetic equations [ jl6| . 
It has been shown that they establish the complete con- 
servation laws. These A-contributions following from the 
nonlocality of the scattering integral read for the energy, 
pressure tensor and density 

v h =\Y.j [(P+?)iA4,-+(fc-g)iA 3j -piA 2j ] , 
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where * = \T^\ 2 2^5{e 1 +e 2 -e 3 -e i )f } h{l~ f 3 - / 4 ) is 
the probability to form a molecule during the delay time 
A*. 

While these correlated parts are present in the numer- 
ical results and can be shown to contribute to the con- 
servation laws we will only discuss the thermodynamical 
properties in terms of quasiparticle quantities to compare 
as close as possible with the mean field or local BUU ex- 
pressions. The discussion of these correlated two - parti- 
cle quantities are devoted to a separate consideration. 



B. Dynamical thermodynamical variables 

We want now to construct the time dependent global 
termodynamical variables. From the distribution func- 
tion /(p, r, t) the local density, current and energy den- 
sities are given by 
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which are computed directly from the numerical solution 
of the kinetic equation in terms of test particles. Please 
note that the above kinetic energy includes the Fermi 
motion. 



1 . Temperature 

The global variables per particle number like kinetic 
energy, Fermi energy and collective energy are obtained 
by spatial integration 
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where we have used the local density approximation 124] . 
Now we adopt the picture of Fermi liquid theory which 
connects the temperature with the kinetic energy as 

£ K {t) = \s F {t) + £ co n(t) + -^^T(t) 2 (15) 
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from which we deduce the global temperature. The def- 
inition of temperature is by no means obvious since it 
is in principle an equilibrium quantity. One has several 
possibilities to define a time dependent equivalent tem- 
perature which should approach the equilibrium value 
when the system approaches equilibrium. In ]25|,^6| the 
definition of slope temperatures has been discussed and 
compared to local space dependent temperature fits of 
the distribution function of matter. This seems to be a 
good measure for higher energetic collisions in the rela- 
tivistic regime. Since we restrict here to collisions in the 
Fermi energy domain and do not want to add coalescence 
models we will not use the slope temperature. Moreover 
we define the global temperature in terms of global en- 
ergies which are obtained by local quantities rather than 
defining a local temperature itself. This has the advan- 
tage that we do not consider local energy fluctuations but 
only a mean evolution of temperature. 



and g the density evolution. We see that depending on 
the considered volume sphere we obtain different global 
densities. We follow here the point of view that the mean 
square radius will be used as a sphere to define the global 
density. This is also supported by the observation that 
the mean square radius follows the visible compression. 
This becomes evident in figure || for a symmetric reac- 
tion at higher energies where at 40fm/c we see a clear 
compression. If we define the volume by a density cut 
n > uq/10 in spatial domain we will not see compression 
at all since the matter is evaporating and this volume 
increases correspondingly to compression. Therefore we 
think that the sphere with the mean square radius is a 
good compromise. 

III. ISO-NOTHING CONDITIONS IN 
EQUILIBRIUM 



2. Energy and pressure 



The mean field part of the energy is given by 
U{t) = £^{t) - S K (t) - £ Bmn (t) 
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from which one deduces the pressure per particle 



(16) 



2n (s + l)ng 



drn(r, t) 



(17) 



In order to compare now the local BUU with the non- 
local BUU scenario we consider the energy which would 
be the total energy in the local BUU without Coulomb 
energy 



£(t)=£ K (t)-£ con (t) + U(t). 



(18) 



This expression does not contain the two - particle cor- 
relation energy which is zero for BUU and the Coulomb 
energy. The reason for considering this energy for dy- 
namical trajectories is that we want to follow the tra- 
jectories in the picture of mean field and usual spinodal 
plots. 



3. Density 

To define the density exhibits to some extend a prob- 
lem. To illustrate this fact we have plotted in figure |l| 



Let us first recall the figures of mean field isotherms 
in equilibrium. The mean field Skyrme and Born corre- 
lational energy is 
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with the kinetic energy in terms of standard Fermi inte- 
grals and the density 
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with g the spin, isospin,... degeneracy. The correspond- 
ing pressure reads 
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We obtain the typical van der Waals curves in figure 
^. Since we have neither isothermal nor isochoric nor 
isobaric conditions in simulations, shortly since we have 
iso-nothing conditions, we have to find a representation 
of the phase transition curves which are independent of 
temperature but which reflects the main features of phase 
transitions. This can be achieved by the product of en- 
ergy and pressure density versus energy density in figure 
H below. This plot shows that all instable isotherms ex- 
hibit a minimum in the left lower quarter. There the en- 
ergy is negative denoting bound state conditions but the 
pressure is already positive which means the system is un- 
stable. The first isotherm above the critical one does not 
touch this quarter but remains in the right upper quarter 
where the energy and pressure are both positive and the 
system is expanding and decomposing unboundly. The 
left upper quarter denotes negative pressure and energy 
indicating that the system is bound and stable. 
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FIG. 1. The time evolution of central collisions Ni + Au at 25MeV per nucleon. The density contours in the spatial plane 
(x,0, Z) is plotted in the left figure in the range — 15, 15fm and the arrows characterize the values of the local current J 
according to (p^[). The middle figure gives the density profile in fm -3 of the beam direction (solid line) and perpendicular to 
the beam (dotted line) . The corresponding right panel shows the global density ratio to nuclear density no = 0.16fm~ 3 defined 
in a sphere versus the radius of the sphere. The mean square radius is marked explicitly by a dot on the radius axes. 
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FIG. 2. The time evolution of central collisions X e + Sn at 50MeV per nucleon analogous to figure 
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FIG. 3. The isotherms for the pressure density versus vol- 
ume (above) and for the product of pressure and energy den- 
sity versus energy density (below). The temperatures are 
T= 1,4, 7, 10,13, 16, 19MeV. 

In order to achieve now a temperature independent 
plot we scale both axes of figure || (below) with a temper- 
ature dependent polynomial and achieve that all critical 
isotherms are collapsing on one curve in the left lower 
quarter, see figure [| 
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FIG. 4. The isotherms of the product of pressure and 
energy density versus energy density scaled by a tem- 
perature polynomial f(T), g(T). The temperatures are 
T = 1,4, 7, 10, 13, 16, 19MeV respectively. All critical 
isotherms collapse on one line in the left lower quarter. 

The first isotherms above the critical one does not en- 
ter the left lower quarter. We consider this scaling as 
adequate for iso-nothing conditions. A phase transition 



should be possible to observe if there occurs a minimum 
in the left half of this plot at negative energies. 

The idea of plotting combinations of pressure and en- 
ergy is similar to the one of softest point [^7| in analyzing 
QCD phase transitions. There the simple pressure over 
energy ratio leads to a temperature independent plot due 
to ultrarelativistic energy - temperature relations. In our 
case we have a Fermi liquid behavior at low temperatures 
and have to scale differently. In particular we have used 
in figure |i| the temperature dependent polynomials 

g[x = T/MeV] = 1.2 + -lQ~ l x + -10 -3 a: a + -KrV 
3 3 5 

+ -10~ 6 x 6 
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f[x = T/MeV] = 1.2 + -lQ~ 2 x + -KTV + 10 _5 a: 4 
3 3 

H-KTV (22) 

which are producing a temperature independent plot in 
figure |] for the specific used mean field potential param- 
eterization. 



IV. NONEQUILIBRIUM THERMODYNAMICS 

Let us now inspect the dynamical trajectories for the 
above defined temperature, density and energy. In figure 
H we have plotted the dynamical trajectories for a charge 
- symmetric reaction of Xe on Sn at 25MeV lab - energy. 
The solution of the nonlocal kinetic equation is compared 
to the local BUU one. One sees in the temperature versus 
density plane that the point of highest compression is 
reached around 60fm/c with a temperature of 9 MeV. 

After this point of highest overlap or fusion phase we 
have an expansion phase where the density and temper- 
ature is decreasing. While the compression phase is de- 
veloping similar for the BUU and for the nonlocal ki- 
netic equation we see now differences in the development. 
First the temperature of the nonlocal kinetic equation is 
around 2MeV higher than the local BUU result. This is 
due to the release of correlation energy into kinetic energy 
which is not present in the local BUU scenario. After this 
expansion stage until times of 120fm/c we see that the 
BUU trajectories come to a rest inside the spinodal re- 
gion while the nonlocal scenario leads to a further decay. 
This can be seen by the continuous decrease of density 
and increase of energy. Since matter is more decomposed 
with the nonlocal kinetic equation we also heat the sys- 
tem more due to Coulomb acceleration. This leads to the 
enhancement of temperature compared to BUU. An os- 
cillating behavior occurs at later times which reflects an 
interplay between short range correlation and long range 
Coulomb repulsion. The decomposition leads almost to 
free gaseous matter after 300fm/c as can be seen in the 
energy versus density plot. 

Please note that although the trajectories seems to 
equilibrate inside the spinodal region when one considers 
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the temperature versus density plane, we see that in the 
corresponding energy versus temperature plane the tra- 
jectories travel already outside the spinodal region. This 
underlines the importance to investigate the region of 
spinodal decomposition in terms of a three dimensional 
plot instead of a two dimensional one like in the recently 
discussed caloric curve plots. Different experimental sit- 
uations lead to different curves as long as the third coor- 
dinate (pressure or density) remains undetermined. 

The iso-nothing plot analog to figure |J in the lower 
left corner shows that the point of highest compression is 
linked to a first instability seen as a pronounced minimum 
of the trajectory in the left quarter. This is connected 
with a pronounced surface emission and connected with 
anomalous velocity profiles j28). We will call this phase 
surface emission instability further on. At 180 fm/c we 
see a second minimum which is taking place inside the 
spinodal region. This instability we might now attribute 
to spinodal decomposition since the trajectories develop- 
ing slower and remain inside the spinodal region. The 
BUU shows the same qualitative minima but the matter 
rebounds and the trajectories move towards negative en- 
ergies again. In opposition the nonlocal scenario leads to 
a further decomposition of matter as described above. 

In the next figure |(] we have plotted the same reac- 
tion as in figure || but at higher energy of 33MeV. We 
recognize a higher compression density and temperature 
than compared to the lower bombarding energy. Conse- 
quently the trajectories develop further towards the un- 
bound region of positive energy after 300fm/c. While the 
first surface emission instability is strongly pronounced 
we see that the second minimum in the iso-nothing plot 
is already weaker indicating that the role of spinodal de- 
composition is diminished. The trajectories in the tem- 
perature versus density plot comes still in the spinodal 
region at rest but travels already outside the spinodal re- 
gion if the energy versus temperature plot is considered. 
This shows that the trajectories start to develop too fast 
to suffer much spinodal decomposition. 

If we now plot the same reaction at 50MeV in figure 
we see that the trajectories come at rest outside the spin- 
odal region whatever plot is used and no second minima 
is seen anymore in the iso-nothing plot. But, the surface 
emission instability is still very pronounced and is proba- 
bly here the leading mechanism of matter disintegration. 

We might now search for a situation where we have the 
opposite extreme that is we search for a reaction with as 
less as possible surface emission instability and as much 
as possible spinodal decomposition. For this reason we 
might think on asymmetric reactions since the different 
sizes of the colliding nuclei might suppress the surface 
emission. Indeed as can be seen in figure || for an asym- 
metric reaction of Ni on Au at 25MeV lab-energy with 
nearly the same total charge as in the reaction before that 
the surface emission instability is less pronounced while 
the spinodal instability is much more important. There 
appears even a third minima showing that the matter 
suffers spinodal decomposition perhaps more than once 



if the bombarding energy is low enough and a long oscil- 
lating piece of matter is developing. 

The higher bombarding energies now show the same 
qualitative effect in that it pronounces the surface emis- 
sion instability and reduces the importance of the spin- 
odal decomposition as can be seen in figures | and 0. 
Please note that much smaller compression densities and 
temperatures are reached in these reactions compared to 
the more symmetric case of Xe on Sn. 

V. SUMMARY 

The nonlocal kinetic theory leads to a different 
nonequilibrium thermodynamics compared to the local 
BUU. We see basically a higher energetic particle spec- 
tra and a higher temperature of 2MeV. This is attributed 
to the conversion of two-particle correlation energy into 
kinetic energy which is of course absent in local BUU 
scenario. 

By constructing a temperature independent combi- 
nation of thermodynamical variables we are able to 
investigate the signals of phase transitions under iso- 
nothing conditions. Two mechanisms of instability 
have been identified: surface emission instability and 
spinodal decomposition. We predict for the cur- 
rently investigated reactions seen in table [j] which ef- 
fect should be the leading one for matter decomposition. 





25 MeV 


33 MeV 


50 MeV 


IlNi + ^| 7 Au 
If Xe + Ih'Sn 


S 

c s 


C S 
C(S) 


C(S) 

c 






15 MeV 


33 MeV 


60 MeV 


JfGd + gfU 
^Ta + « 7 Au 


C S 


C S 
C(S) 


C 
C 



TABLE I. The prediction of the leading mechanisms of 
matter disintegration for two reactions with equal total charge 
but asymmetric entrance channels. Surface compression is de- 
noted by C and spinodal decomposition by 5*. 



In the reactions with bombarding energies higher than 
the Fermi energy the fast surface eruption happens out- 
side the spinodal region. For even higher energies there is 
not enough time for the system to rest at the spinodal re- 
gion. The trajectories simply move through the spinodal 
and the system decays before it comes to an equilibrium 
- like state inside the spinodal region. 
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